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ABSTRACT 

The Local Group is often seen to be a quiescent environment without significant merger 
events. However an ancient major merger may have occurred in the most massive 
galaxy as suggested by the M31 classical bulge and its halo haunted by numerous stellar 
streams. Numerical simulations have shown that tidal tails formed during gas-rich 
major mergers are long-lived and could be responsible for old stellar streams and likely 
induce the formation of tidal dwarf galaxies (TDGs). Using several hydrodynamical 
simulations we have investigated the most prominent tidal tail formed during the first 
passage, which is gas-rich and contains old and metal poor stars. We discovered several 
striking coincidences after comparing its location and motion to those of the Milky 
Way (MW) and of the Magellanic Clouds (MCs). First, the tidal tail is sweeping a 
relatively small volume in which the MW precisely lies. Because the geometry of the 
merger is somehow fixed by the anisotropic properties of the Giant Stream (GS), we 
evaluate the chance of the MW to be at such a rendez-vous with this gigantic tidal tail 
to be 5%. Second, the velocity of the tidal tail matches the LMC proper motion, and 
reproduce quite well the geometrical and angular momentum properties of the MW 
dwarfs, i.e. the so-called disk of satellites, better called Vast Polar Structure (VPOS). 
Third, the simulation of the tidal tail reveals one of the formed TDG with mass and 
location almost comparable to those of the LMC. Our present modeling is however 
too limited to study the detailed interaction of gas-rich TDGs with the potential of 
the MW, and a complementary study is required to test whether the dwarf intrinsic 
properties can be accounted for by our scenario. Nevertheless this study suggests a 
causal link between an expected event, an ancient, gas-rich major merger at the M31 
location, and several enigma in the Local Group, namely the GS in the M31 outskirts, 
the VPOS almost perpendicular to the MW disk, and the presence of the MCs, two 
Irr galaxies near the MW. 
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1 INTRODUCTION 

Despite their proximity, the properties of MW dwarf galax- 
ies are far from being unanimously interpreted. They are 
often considered as remnants of primordial galaxies that 
would have escaped to the numerous merger events expected 
in the hierarchical scenario. Within the ACDM paradigm 
they are believed to be dark matter dominated sub haloes 
that are gravitationally trapped in the MW halo and have 
not y et merged completely with their host l|White fc Reesl 
Il978l h As the ACDM predicted a considerable number of 
massive sub-haloes relatively to the small number of ob- 
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served dwarf galaxies surrounding the MW and other giant 
spirals, this leads to objects considerably dominated by their 
dark matter content. Such objects, including UFDs, are cur- 
rently followed to search for dark matter annihilation, with- 
out any positive detection yet (see e.g., Ackermann et all 



|2012| ; iGeringer-Sameth fc Koushiappadl201lD . 

However some dwarf galaxies may be tidally formed 
during the numerous m erger events ex pected in the frame 
of the ACDM scenario |Kro upa 2012) , leading to o bjects 
mainly free of dark matter. lOkazaki fc Taniguchil (|2000h 
claimed that most dwarf galaxi es could b e of t idal ori- 
gin, although thi s is disputed by iBournaudl l|20ld ) and by 
IWen et all l |2012t ). A full description of the number density 
of TDGs is still lacking, as it has to account for gas-rich 
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mergers that likely occurred in the past, as well as to con- 
sider a representative sampling of (major) merger orbital 
parameters and to describe the production and destruction 
rate with time of TDGs within evolving tidal tails. 

A significant part of our knowledge of dwarf galaxy 
properties is based on the Local Group content, includ- 
ing irregular dwarf galaxies (dlrr) with the two Magellanic 
Clouds (MCs), spheroidal dwarf galaxies (dSph), and Ultra- 
FaintDw arf galaxies (UFD), which have been recently dis- 
covered JWillman et al. 2005; Sak amoto fc Hasegawa 200r3; 



Zucker et al.ll2006l: iBelokurov et al.ll2006l . l2007l : IWalsh et al l 
20071 ). I Einasto et al. I (|l974r i pointed out that dSphs lie close 



to the MW (mean distance 193 kpc), whereas the dlrr are 
more distant (572 kpc), except the MCs. In o rder to e xplain 
this distance-dependent morphological bias, iMaveij (|201lh 
suggested that the closest dwarf galaxies undergo physical 
phenomenons that changed their morphology from dlrr to 
dSp h, such as stripping, stir ring or tidal force. 

iKallivavalil et alJ l|2009h have re-estimated the proper 
motion of the Large Magellanic Cloud (LMC) to a signifi- 
cantly larger value than previous estimates, i.e. a velocity of 
378 km/s relatively to the MW. Th us either MCs are passing 
for the first time close to the MW l|Besla et al.ll2007f ) or one 
needs to assume a quite large dark matter content for the 
MW with M 30 ofc pc = 2.7 ± 0.5 x 10 12 M (jWatkins et all 
2010), i.e. a value 50 times larger than its baryonic content 
and even larger to tha t of M31. A recent measurement of the 
LeoII proper motion ijLepine et al provides a similar 

constraint, letting open the question of whether some of the 
MW dwarf galaxies may be unbound, and as such, would 
not be MW satellites. 

In addition, most of the MW dwarf galaxies seems to be- 
long to a plan-like structure, the Plane of Satellite (VPOS) 
which is found to be roughly perpendicular t o the MW disk 
Kunkel fc Demerdll97lSULvnden-Bel]|ll97rJ ). iKroupa et all 



2005h updated this idea by using the 11 classical dwarf 



galaxies and compared the satellite spatial distribution with 
the isotropic spatial distribution expected from ACDM sim- 
ulations. They concluded that the dwarf galaxy distribution 
cannot derive from the expected ACDM distribution, and 
iMetz et al] l)2007ah deduced that their spatial distribution 
is neither spherical nor mildly prolate with a confidence of 
more than 99%. Moreover, the VPOS could be a somewhat 
permanent structure as the angular momenta of five dSphs 
plus the MCs ind icate a coherent motion within the VPOS 
ilMetz et al 1l2008t^_ 



Met z et all (Hb09b) summarized the different possible 



concepts to explain the concordance between the spatial and 
angular momentum distribution of the VPOS galaxies. They 
investigated the scenario of a dwarf g alaxy group surround- 
ing the M Cs (iD'Onghia fc Lakell2008I). possibly comi ng from 
a filament l)Zentner et al .1120051 : lLibeskind et al J2005h , which 
would be entering the MW halo in ordered motions. None 
of these explanations seem fully satisfactory, because it re- 
quires a quite unexpected compact group as a progenitor, 
and because the s uper galactic plan e is almost perpendicu- 
lar to the VPOS (|Metz et alj|2009bl ). 

A fully different alternative would be that MW dwarf 
galaxies originate from an ancient tidal tail (see a detailed 
discussion on the possible tidal origin in iKroupal 1 19971 ; 
IMetz et al. 2009b), which would explain their ordered mo- 
tion. In fact such an alternative was suggested quite early 



l|Lvnden-Belll Il976h and iPawlowski et al] l|201lf ) have re- 
cently investigated the possibility of a major merger or a 
fly-by interaction (mass ratio 1:1 or 4:1) that would have oc- 
curred in the early history of the MW, producing tidal tails 
and then TDGs. Assuming an orbital angular momentum 
aligned to the VPOS, the tidal tail associated to this early 
interaction in the MW history would have formed TDGs 
that could further populate the VPOS. While this is an in- 
teresting suggestion, perhaps the MW is not the best candi- 
date for being a merger remnant, w hich are generally associ- 
ated with classical bulge galaxies (Kormcndy & Kcnnicutt 
2004). Moreover, as mentioned by IPawlowski et al] (|20 1 IT ) . 
the fly-by alternative would have let the interloper well de- 
tectable, at the fringes of the Local GroupQ. 

All the above ideas are very imaginative and they may 
resolve the origin of the VPOS and why dwarf galaxies are 
part of it. However this is at the cost of introducing an- 
other ad-hoc assumption. Alternatively, we may consider the 
past history of the MW in the whole context of the Local 
Group, i n which the bar y onic c ontent is dominated by M31. 
Quoting Ivan den Berghl (|2005h : 1 'Both the high metallicity 
of the M31 halo, and the r 1//4 luminosity profile of the An- 
dromeda galaxy, suggest that this object might have formed 
from the early merger and subsequent violent relaxation, 
of two (or more) relatively m assive metal-rich ancestra l ob- 
jects". As further noticed by lYang fc Hammer! (|2010h . the 
disk plane of M31 is seen almost perpendicular, implying 
that tidal tails possibly formed after a gas-rich merger in 
the M31 history could be part of a hyper-plane that includes 
the MW. By trac i ng ba ck in time the position of the LMC, 
lYang fc Hammed (2010) have shown that the LMC could be 
near M31, 5 to 8 Gyr ago. Moreover. lHammer et all (|2010h 
have proposed that M31 could be the result of a major, 
gas-rich merger because it provides a simple and common 
interpretation of most of its exceptional properties, includ- 
ing the Giant Stream, the outer thick disk, as well as the 
giant ring. Constraints from stellar population ages in these 
different halo substructures imply a first clos e passage 8 to 9 
Gyr ago and a fusion time 5.5 to 6 Gyr ago l|Hammer et al] 
l2010h . 

This paper intends to verify whether or not the ordered 
motion of the MW dwarf galaxies could be entirely due to 
their origin as TDGs formed from a tidal tail caused by an 
ancient interaction in the history of M31. The goal is to 
verify whether this alternative may explain two exceptional 
features in the Local Group: the VPOS and the MW-LMC- 
SMC proximity. iRobotham et"al] ()2012h estimate the occur- 
rence of the MW-LMC-SMC configuration for 414 MW-like 
galaxies m a local volume of 1.8 x 10 5 Mpc 3 (0.01 < z < 
0.055). The chance to found out a galaxy having two close 
companions at least as massive than the SMC within a pro- 
jected separation of 70 kpc and a radial separation of 400 
km.s -1 is only 0.4 % (2 examples). Moreover, in each of these 
configurations there is another L* luminosity galaxy within 
1 Mpc from the MW-like galaxy. This may support the im- 



1 We do not consider the possibility that the interloper could be 
indeed the LMC, because this would require a very problematic, 
huge difference of cosmological growth between LMC and the 
MW. 
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portance of M31 for the MCs formation, and then, account- 
ing for the VPOS, for all the MW dwarves. 

In Sect. 2, we present the VPOS properties and refine 
its statistical significance by including with additional con- 
straints on the dwarf angular momenta. In Sect. 3, we de- 
scribe both the physical and numerical models that we adopt 
for reproducing the VPOS. In Sect. 4, we try to reproduce 
the VPOS with a simple tidal tail toy model. Then in Sect. 
5, the model is i mpro ved to closely match the work made by 
lHammer et al.l (|2010l ). providing more realistic initial con- 
ditions for the tidal tail. Sect. 6 summarizes the observa- 
tional coincidences that support the M31 scenario, discusses 
the possible falsifications, and describes the required further 
studies for testing them. 

In this paper, we use two coordinate systems. One is 
the Galactocentric coordin ate which is centered on the MW 
|van der Marel et al.ll2002h . This system is used in Sect. 2 
and Sect. 4. In Sect. 5, we use a "projection" coordinate 
which is centered on the current M31 position with the z- 
axis pointing from M31 to the Sun, the x- and t/-axis parallel 
to the East and the North of celestial coordinate at M31 
position. In this projection coordinate, MW is located at (x, 
y, z) = (6.42, -2.72, 788.85) kpc and the Sun at (0, 0, 785) 
kpc, assuming that the distance of M31 is 785 kpc. 



2 VPOS PROPERTIES 
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Figure 1. Distribution function of the satcllitc-MW distances. 
Large red points represent the dSphs, small red ones the UFDs, 
while blue squares represent the dlrrs. A line separates the MW 
companions considered in this study (closer than 300 kpc) from 
the other dwarf galaxies. 



Table 1. Basic properties of the 11 classical dwarf galaxies. From 
left to right: galactic longitude /, galactic latitude b (both in de- 
grees), heliocentric distance R (kpc), V-band ab solute mag nitude 
M v , a nd refe rences: (1 i = iMcConnachid <2012h . (2) = iMateol 
l|l998h . (3) = Ivan den Berghl dl99Sh 



2.1 Defining the VPOS dwarf galaxies 

The distances between the dwarf galaxies and the MW can 
be used as a simple criterion to verify whether dwarf galaxies 
are MW companions. The gravitational force of the Galaxy 
becomes weaker for distances larger than the v irial radius, 
which ranges from 150 to 300 kpc. Following iMetz et al.l 
(2009a), we consider as companions only the galaxies in a 
sphere with a radius of 300 kpc (see Fig.Q]). With this simple 
criterion, all the dlrr but the MCs are discarded, while most 
of the dSph and UFDs are taken into account. The next 
goal is to get a complete sample of dwarf galaxies within 
this sphere with a cut in absolute V-band magnitude. Two 
factors make this task difficult: the 11 classical dwarf galax- 
ies were discovered by different sets of observations, and the 
Zone of Avoidance hides a part of the sky. Thus in principle, 
our sample may be affected by strong observational biases. 
The fact that the newly discovered satellites by the SDSS 
are systematically fainter than the classical dwarf galaxies is 
howeve r reassuring. Even Canes Vanetaci I with a Mv= -7.5 
± 0.5 (|Zucker e t al. 2006) is fainter than the faintest clas- 
sical dwarf, suggesting that no classical dwarf could have 
been misse d in the sky not sam pled by the SDSS (80% of 
the sky, see lKoposov et al"1l2008l ). except maybe in the Zone 
of Avoidance. Accounting for the UFDs found by the SDSS 
that includes a significant part of the VPOS could therefore 
bias the resulting statistics. 

Thus, we will consider only the 11 classical dwarf galax- 
ies grouped in a sample named "dSphMC": LMC, SMC, 
Sagittarius, Fornax, Leo I, Leo I I, Sculptor, S e xtans , Ca- 
rina, Draco , Ursa M inor, following lKroupa et al.l ()2005l L and 
IMetz et all (|2007al ). Their basic properties are listed in Ta- 
ble [1] 



Name 


I 


b 


R 


M v 


Ref. 


LMC 


280.46 


-32.88 


50 ± 2 


-18.1 ± 0.1 


(1) 


(3) 


SMC 


302.80 


-44.32 


60 ± 4 


-16.2 ± 0.2 


(1) 


(3) 


Fornax 


237.10 


-65.65 


138 ± 12 


-13.4 ± 0.3 


(1) 


(2) 


Leol 


225.98 


49.11 


250 ± 15 


-12.0 ± 0.3 


(1) 


(2) 


LeoII 


220.16 


67.22 


210 ± 14 


-9.8 ± 0.3 


(1) 


(2) 


Sextans 


243.49 


42.27 


86 ± 5 


-9.3 ± 0.5 


(1) 


(2) 


Carina 


260.11 


-22.22 


102 ± 6 


-9.1 ± 0.5 


(1) 


(2) 


Draco 


86.36 


34.72 


82 ± 6 


-8.8 ± 0.3 


(1) 


(2) 


UMi 


104.95 


44.80 


66 ± 7 


-8.8 ± 0.5 


(1) 


(2) 


Sgr 


5.60 


-14.08 


24 ± 2 


-13.5 ± 0.3 


(1) 


(2) 


Sculptor 


287.53 


-83.15 


87 ± 5 


-11.1 ± 0.5 


(1) 


(2) 



2.2 Spatial properties of the VPOS 

We follow the methodology of IMetz et all l|2007al . l2009aT ) 
who have studied in detail the spatial properties of the 
VPOS, by fitting the position of the VPOS using a least- 
square method. They investigated the uncertainties associ- 
ated to the direction perpendicular to the VPOS due both to 
the uncertainties on the distances between the MW and the 
satellites (using Monte-Carlo simulations) and to the fitting 
model (using bootstrap resampling). 

Using the data listed in Table Q] we confirm that the 
dominant source of uncertainty is due to the fitting model 
(see Table [2j|. Figure [2] shows the resulting VPO S fit, which 
is foun d to be very similar to that derived by IMetz et al.l 
l|2007ah . 
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Table 2. Results of the VPOS fit. From left to right: galactic 
longitude and latitude of the vector perpendicular to the VPOS, 
lypos an d °VPOS (deg), uncertainty due to the plane fitting 
method, Errorl (deg), uncertainty due to the uncertainty associ- 
ated to the VPOS dwarf galaxy distance, Error 2 (deg), minimal 
distance from the MW center to the VPOS plan, D-center (kpc), 
standard deviation of the minimal distance from the VPOS dwarf 
galaxy position to the VPOS, thickness, T (kpc). 



Wpos 


WpOS 


Errorl 


Error2 


D-Center 


T 


157.4 


-12.5 


16.4 


1.03 


8.07 


18.5 



2.3 Kinematic properties of the VPOS 

In addition to spatial correlations, iMetz et al.l (|200St ) have 
robustly confirmed that most of the individual dwarf mo- 
tions lie close to the VPOS. We have extended and con- 
firmed their study by adding two dwarf galaxies, i.e., LeoII 
and Sextans (see Table[3]). In the left panel of Fig. we show 
that the angular momentum directions of all the dwarf ex- 
cept Ursa Minor and Sagittarius, lie close to the uncertainty 
range of the VPOS perpendicular direction (notice that Ca- 
rina and Draco are only marginally consistent). Besides this, 
the Sculptor motion is found to be lying in the VPOS but 
counter-rotating. Indeed, its angular momentum is shifted 
by 180°. 

As the MW dwarf galaxies orbit within a thick plane 
that is not rotatio nally supported, the t erm Vast Polar 
Structure (VPOS) (|Pawlowski et all l2012h is better than 
Disk of Satellite to describe the true nature of this struc- 
ture. 



2.4 Significance of the VPOS 

We test the hypothesis of an isotropic distribution for the 
10 classical dwarf galaxy positions and velocities (Leo I is 
discarded because its velocity is unkn own). To do this, w e 
used the statistical test proposed by IMetz et all (l2007al ). 
which is based on fitted planes to isotropic distributions (see 
details in Appendix). The null hypothesis is rejected with a 
probability of 99.3 % and 97.0 % for the positions and the 
velocities, respectively. In an isotropic model, the position 
and velocity distributions are independant, which leads us to 
reject the isotropic distribution hypothesis in the full phase 
space with a probability of 99.98 %, which further supports 
the significance of the VPOS. 

However more realistic cosmological simulations are 
providing more compl ex satellite geometries around their 
host galaxy (see e . g. | Libeskind et all 120051 ; iDeason et all 



l2012h . Deason et all l|2012l) have used simulations that avoid 



an over-efficient tidal stripping of the satellites before fusion, 
generating a signi ficant fraction of reconstructed spirals af- 
ter major fusions (|Font et al.ll201lf) . These simulations are 
thus more realistic to figure out what is the distribution of 
satellites in host galaxies. Up to 20% of the hosts may show 
a polar distribution of their satellites, while less than 10% 
of the satellites may have some coherent motions. The MW 
satellites being both pol ar distributed and w ith coherent mo- 
tions within the VPOs TDeason et all |2012r ) concluded that 
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Figure 2. VPOS position derived from the 11 classical dwarf 
galaxies in the MW outskirt. Red points represent the dSphs 
galaxies, while blue squares represent the MCs. The red lines rep- 
resent the uncertainties on the heliocentric distance of the dwarf 
galaxies. The MW is indicated by a 100 kpc line. Top panel: edge- 
on view of the VPOS (magenta dashed line). Bottom panel: same 
as the top panel but rotated by 90° along the polar-axis of the 
MW disk. 



"if a substantial number of the classical dwarf galaxies of the 
Milky Way (e.g. 7 out of 10) were found to have orbital poles 
aligned with the normal to the disc of satellites, then this 
would be inc onsistent wi t h the results of our simulations" . 
As shown bv lMetz et"al] l|2007al ) this is because the VPOS 
is particularly thin which requires the accretion of a large, 
hypothetic group of dwarf galaxies that is otherwise natu- 
rally brought by remnant tidal tails. Thus additional data 
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on the proper motions are still required, though the present 
data favor an exceptional MW satellites alignment of their 
locations and motions (see Fig. [2] [3}. 



3 FORMING THE VPOS AS A REMNANT OF 
A TIDAL TAIL 

3.1 Physical model 

lHammer et"all (|2010h have proposed that a gas-rich major 
merger occurred during the history of M31. The y used an 
N-bo dy, hydrodynamical simulation (GADGET2, ISpringefl 
120051 ) to reproduce most of the M31 structures such as the 
bulge, the thin and thick disks, the Giant Stream, and the 
10 kpc ring. In their model, the first passage between the 
galaxies occurred more than 8 Gyr ago, and the fusion time 
2.5-3 Gyr later. A significant part of the matter ejected dur- 
ing the first-passage lies in a tidal tail. As we will show in 
Sect. 15. H for most of the models reproducing M31, this tidal 
tail is found to pass close to the MW. 

Several observations and simulations have suggested 
that dwarf galaxies, named Tid al Dwarf Galaxies (TDGs) , 
can form in such tidal t ails (seelOkazaki fc Taniguchill2000l . 
and references therein). iBournaud fc Dud (|2006h have run 
96 major merger simulations in order to better under- 
stand TDG properties and their formation mechanisms. 
Most of their simulations revealed the formation of TDGs 
with stellar masses larger than 10 8 Mq. However, amongst 
the MW companions, 9 of the 11 classical dwarf galax- 
ies have stellar mas ses ranging from ~ 10 s to ~ 10 7 Mq 
JStrigari et al. 20081) . Using simulations with higher resolu- 
tion. IBournaud et alj (|2008l ) have further shown that young 
stellar objects, star clusters, and TDGs with masses rang- 
ing from 10 5 to 10 9 M can also be reproduced. All these 
simulations suggest that it is possible to form TDGs with a 
relatively large range of mass, as a result of a major merger 
event. 

In the following sections, we investigate whether the 
remnant of such a ti dal tail could be at the origin of the 
VPOS dwarf galaxies. IPawlowski et all l|201ll ) were the first 
to study in detail such a hypothesis. They simulated the 
evolution of tidal tails created by major mergers or fly-byes 
in the MW past, ~ 10 Gyr ago. They explored the orbital 
parameters of the progenitors to match the properties of 
the VPOS. However, a fly-by seems unlikely, since no other 
massive galaxy but M31 is observed in the Local Group. 
In the present study, we investigate the formation of the 
VPOS as a result of a tidal tail ejected during a merger that 
occurred in the p ast history of M3 1 rath er than in that of 
the MW. Indeed, lYang fc Hammerl (|201Gu have shown that 
the position of the LMC could be traced back toward M31, 
~ 6 Gyr ago, assuming an appropriate choice for the M31 
tangential velocity and the MW halo shape. This suggests 
that all VPOS dwarf galaxies could have formed as TDGs 
resulting of a major merger event in the M31 history. 



3.2 Numerical models 

To investigate this scenario, we developed a C code that sim- 
ulates the trajectories of tidal tails populated with TDGs. 



They are formed and ejected during a major merger, and 
eventually interact with the MW. 

The small mass of the TDGs, as well as the large 
distance between M31 and the MW, make the tidal ef- 
fects undergone by the two spiral galaxies negligible. 
This allows us to assume that their gravitational poten- 
tials are static, and to model them analytically. Following 
I Yang fc Hammed j2010h. we assume Navarro-Frenk- White 
halos jNavarro et alj 19971 ') , Hernquist bulges jHernquistl 
1990), and exponential disks ijYang fc Hammerl l2010l ) for 
M31 and the MW, while dwarf galaxies and the LMC are 
assumed to be point masses. 

In the simulations, we first trace back the position of 
the LMC, the MW, and M31 backwards in time until the 
LMC reaches the outskirt of M31. The unknown tangential 
velocity of M31 is taken as a free parameter, while other 
parameters are listed in Table [T]& Table [3] Since the LMC 
is assumed to be a TDG remnant, its mass is just equal to 
its baryonic mass, ~ 3.10 9 Mq. We have verified that the 
dynamical friction is negligible in the computation of the 
LMC trajectory. 

Second, we used numerical models of tidal tail formation 
to set up initial conditions for the positions and velocities 
of the TDGs in the surrounding of M31. The position of 
the LMC back in time in the surrounding of M31 is used to 
constrain the tidal tail properties. Finally, we let evolve the 
TDGs until present-time. At the end of the simulation, we 
analyze the TDG positions and angular momenta, and test 
whether they match the present-day VPOS properties. 



4 A TIDAL TAIL TOY MODEL 
CONSTRAINED BY THE LMC 
TRAJECTORY 

4.1 Initial conditions for the TDGs 

The goal of the proposed toy model is to verify whether or 
not the gross properties of the VPOS can be described by 
particles linked to the over-densities found in tidal tails. As 
described in Sect. 13.21 we first trace back in time the position 
of the LMC, MW, and M31, 5.5 Gyr ago. This corresponds 
to the epoch when the LMC was found to be at 50 kpc from 
M31 according to lYang fc Hammerl [|2010h (see their Table 
3). To be consistent with their study, the MW and M31 total 
masses are respectively choosen to be 10 12 Mq and 1.6 x 10 12 
M , while the LMC veloci ty is taken from the study made 
bv lKallivavalil et all |2009l '). i.e., V = 384 km/s. 

The position and velo city profiles of the tida l tail are 
extracted from the work of IWetzstein et all l|2007t ). Indeed, 
they studied the formation of TDGs produced in a major 
merger and estimated the effects of initial conditions on the 
TDG properties. Two snapshots separated in time by 262 
Myr are extracted from their simulation EG2 (see Fig. S}. 
Five density peaks are found to be prominent within the 
tidal tail. We assume that the third peak represents the 
LMC location in the past. This arbitrary choice implies that 
the part of the tidal tail which is located farther than the 
LMC position will escape from M31, which is quite expected 
during the typical evolution of a tidal tail. 

We use the relative positions of the five prominent peaks 
to set up initial positions of five groups of TDGs. Each 
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Table 3. Proper motions of 10 classical dwarf galaxies. From left to right: Radial velocity (km/s), proper motion (mas/yr), both in the 
h eliocentric res t frame , total velocity, radia l veloc ity, tra nsverse velocity (all in km/s in the galacto - centr ic restfr ame), and ref e rence s: (1) 
= Vieira et all d2010h. ( 2) = |Piatek et"al] d2007h. ( 3) = lLepine etHl d201lfl (4) = IWalker etaD ll2008l) . (5) = IPiatek et~aH (|2003h . (6) 



Scholz fc Irwinl dl994f) . (7) = |Piatek et alj <200g) . (8) = llbata et all dl9^7H . (9) = |Piatek et alj | |2006h . 



Name 


V ra d 


V>6 


fi a cos(5) 


V tot 


Vrad 


Vtrans 


Ref. 



LMC 278 0.39 ±0.27 -1.89 ±0.27 347±gg 103^" 331±gg (1) 

SMC 158 -1.01 ±0.29 0.98 ±0.30 221+\l° 29 -i4 219 -i2i W 

Fornax 53 -0.36 ±0.041 0.476±0.046 188±i^ -31±l 185±H (2) 

LeoII 79 -0.033±0.151 0.104±0.113 2681^ 22+^ 267±^ (3) 

Sextans 224 0.10 ±0.44 -0.26 ±0.41 260^^2 88^30 241 -llo (4) 

Carina 229 0.15 ±0.09 0.22 ±0.09 131+^ 27+^ 128+j^ (5) 

Draco -292 1.1 ±0.3 0.6 ±0.4 563^^0 -62l^ 559^4 (6) 

UMi -247 0.22 ±0.16 -0.50 ±0.17 Wl+H -76±t 141 ±58 ( 7 ) 

HO 100 + I ncn+ll 



Sgr 140 -0.88 ±0.08 -2.65 ±0.08 303+^ 138+^ 269+^ (8) 

H69 
-72 



Sculptor 110 0.02 ±0.13 0.09 ±0.13 227±H 79t% 212+^1 (9) 




6 




Figure 3. Left panel: Aitoff map showing the direction of the angular momenta of the classical dwarf galaxies (squares), and their 
uncertainties (lines). The solid ellipse shows the uncertainties associated to the direction perpendicular to the VPOS. The dashed ellipse 
represents the uncertainty associated to the direction of the mean angular momentum. Right panel: (7, Q density map generated using 
Monte-Carlo simulations of isotropic satellite distributions (see Appendix A). A large value of 7 and £, such as found for the MW 
position and velocity distributions, means that the vector perpendicular to the fitted plane have a clear preferred direction and that its 
uncertainty around this direction is small, which means that the fit is significant. 



group, even that with the LMC, is assumed to contain six 
TDGs, which are modeled as point masses with a total mass 
equals to 10 8 Mq, i.e., close to that of the SMC. We choose 
six particles per peak in order to populate enough the out- 
skirts of the MW at the end of the simulation. The mass 
distribution in each group of particles is assumed to follow 
a Plummer density mass profile with a typical radius of 10 
kpc, which constrains the internal velocity dispersion to be 
~ 16 km/s in order to have a static group. The bulk velocity 
of the TDG groups is roughly estimated by comparing the 
positions of the peaks in the two snapshots (see Fig. [4}. The 
tidal tail plane that encompasses the positions and motions 



of the TDGs is rescaled, rotated, and translated in order 
for the third density peak to match both the position and 
velocity of the LMC, 5.5 Gyr ago. 



4.2 The results 

A set of 100 simulations has been run by randomly generat- 
ing internal position and velocity for the TDGs in the five 
groups, consistent with a Plummer model. For each simu- 
lation, we only kept particles for which the trajectories end 
within 300 kpc from the MW centre, i.e., the distance limit 
used in Sect. 12. ll to define a dwarf galaxy as a MW compan- 
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Figure 4. Trajectories of five over-dens ities (black lines) in a 
tidal tail during a time step of 262 Myr llWetzstein et al.|[2007t) . 
The third over-density is assumed to represent the position of 
the LMC back in time. The distance between the LMC and the 
merger remnant center is 50 kpc. 



ion. The VPOS properties and the mean angular momentum 
direction are calculated following the same method that the 
one used for the observations (see Sect. 12.2]) . The resulting 
plane thickness distribution for the 100 simulations has an 
average of 18. 3± 6 kpc, which is consistent with the observed 
value (see Table 

Figure [5] shows typical trajectories of simulated TDGs. 
Particles are found to travel on both sides of the MW disk 
plane, although most of them, including the LMC, pass be- 
low it. There are some particles passing above the MW disk 
plane, with an angular momentum direction shifted by 180° 
relative to the simulated VPOS direction. This configuration 
is similar to the counter-rotation of Sculptor in the MW en- 
vironment. 

In Fig. [6] we plot for the 100 simulations the resulting 
mean angular momentum directions and the directions per- 
pendicular to the simulated VPOS. These directions pop- 
ulate the region occupied by the uncertainty range corre- 
sponding to the observed dwarf galaxies. We find that the 
mean angular momentum direction of the 100 simulations 
is consistent within uncertainties with the mean direction of 
the observed angular momentum. The mean perpendicular 
direction for the 100 simulated VPOS is also close to that 
derived from observations, although not strictly within the 
uncertainty range. However, for some simulations, the per- 
pendicular direction is encouragingly found to lie within this 
range. 

This very simple model illustrates how TDGs associ- 
ated with a tidal tail can easily reproduce most of the VPOS 
properties. By varying the radius of the Plummer model and 
the LMC position in the tidal tail, one can reproduce the 
fraction of dwarf galaxies having the same angular momen- 
tum direction that the LMC (8 dwarf galaxies) and those 
having it inverted as Sculptor, i.e. a very similar optimiza- 
tion that done bv lPawlowski et al.l (|201ll ). 

It is clear that such a model is quite ad-hoc since the 
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Figure 5. Comparison between the positions of the simulated 
TDGs (diamonds) and the observed satellites (squares). The up- 
per and lower panels show the projection in Y-Z and X-Y galactic 
coordinate plane, respectively. The MW is represented as a red 
spiral. The grey lines represent the trajectories of each simulated 
TDG. The in-set panels show a global view of the trajectories of 
the simulated TDGs (grey lines) as well as that of M31 (black 
line). The trajectory of the LMC is also shown as a red line. The 
green line illustrates a counter-rotating trajectory, which may ex- 
plain the motion of Sculptor. 



dominant motion of the simulated TDGs is that of the LMC, 
which is a part of the VPOS. The LM C velocity measure - 
ment used is larger than that found bv lVieira et al.l (|2010h ■ 
Using their value would increase the trace back time needed 
by the LMC to reach M31 by more than 2.5 Gyr. This time 
could match the format ion epoch of the first tidal tail ac- 
cording to the model of lHammer et al l (|2010T ) described in 
Sect.EU i.e 8.5 Gyr ago. Moreover, the TDGs of the VPOS 
remain too close to the LMC and thus do not populate 
enough the region above the MW in Fig. [5] This shows the 
limit of this "toy model" . Finally, such a toy model, as well 
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Figure 6. Aitoff map comparing the VPOS and angular momentum (AM) directions between simulations and observations. The red 
square and the blue dots indicate the direction of the mean AM and the direction of the VPOS, respectively, as derived from observations 
dMetz et al.ll2007al , |200S|'| . Their uncertainties are represented by the red circle (AM) and by the blue dashed circle (VPOS). The open 
triangles and circles represent the directions of the VPOS and the mean AM, as resulting from 100 simulations (see Sect. 14.211 . Their 
mean values are indicated by a green square (VPOS) and a green dot (AM). 



as that proposed by IPawlowski et al.l l|201ll ) , do not allow 
us to distinguish the origin of the tidal tail (i.e, the MW or 
M31). These remarks lead us to make a more realistic model 
of the tidal tail. In the next section, we propose to link far 
more robustly the modeling of M31 and its halo with sim- 
ulations testing the VPOS formation. In principle it should 
bring a number of constraints large enough to rule out or 
support the M31 tidal tail hypothesis. 



5 A TIDAL TAIL GENERATED BY AN 
ANCIENT MERGER AT THE M31 
LOCATION 

5.1 M31 modelling as a remnant of a major 
merger 

lHammer etafl (|20ld ) have run a series of gas-rich, major 
merger simulations designed to reproduce both the M31 
galactic substructures and the extended structures in its out- 
skirts. The orbital geometry of the interaction is close to be 
polar, driven by the need to reproduce the giant 10 kpc ring, 
while stellar ages in the M31 halo substructures constrain 
the pericenter radius. In this family of models, two tidal 
tails systematically form during the first passage, as well 
as additional ones during the second passage and/or fusion. 
Each tidal tail is made by material stripped from each M31 
progenitor, resulting in an angular momentum that com- 
bines that of the corresponding galaxies with the orbital 
angular momentum, which is the dominant component. We 
have t ested a significant series of models from lHammer et al.l 
(2010) (see their Table 3) and found that the first-passage 
tidal tail, associated to the minor encounter, passes system- 



atically close to the MW, i.e., generally at 50 kpc and al- 
ways at less than 200 kpc (see Fig. [8]). In other words, the 
MW lies within the past M31 merger orbital plane, and the 
large area swept by the tidal tail motion, leads to a sig- 
nificant chance for a tidal tail-MW interaction. In fact the 
first-passage tidal tail is the most extended one. It contains 
the largest amount of baryons expelled during the merger, 
and the largest amount of gas, and the lowe st metal abun- 
dance compared to other tidal tails. In the I Hammer et al.l 
(2010) model, the Giant Stream is linked to the tidal tail 
formed during the second passage, and caused by stars re- 
turning to the remnant galaxy leading to the formation of 
loops (see their Fig. 8 as well as a comp lete description of 
the loop mechanism in lWang et alj|2012h . 

Th e model used i n the p resent study is a refined ver- 
sion of lHammer et all (|2010h (10th column in their Table 
3) with a slightly higher mass ratio (3.5 instead of 3), and 
an improvement of the initial conditions, i.e., the two inter- 
lopers at t — are separated by 200 kpc instead of 80 kpc. 
The most significant change is due to the use of a larger 
number of particles (2.4 millions instead of 540 thousands), 
which allows us to generate gravitational instabilities within 
the tidal tail, under the assumption of a softening length of 
0.1 and 0.12 kpc for baryons and dark matter, respectively. 
A more subtle change is about the Giant Stream model- 
ing, which is associated to the s econd order loop inst ead of 
the first order loop (see Fig. 8 in lHammer et alT BoiO). This 
results in a rotation within the orbital plane by approxi- 
mately 50 degrees, which only slightly impacts the relative 
angle between the first tidal tail and the assumed Giant 
Stream. This updated model now reproduces both the Giant 
Strea m and the No r thern loop discovered by the PAndAS 
team iMackev et all l|2010l) , those being related to material 
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stripped in the loop formed at the fusion and at the first 
passage, respectively (see details in Yang et al., 2012, in 
prep.). 

5.2 Initial conditions for the TDGs 



Instea d of relying on the simulatio ns of | Wetzstein et al.l 
|2007l ) as in Sect. [H we now use the lHammer et al.l (|20ld ) 
merger model for M31 (as updated above) to setup initial 
conditions for the position and velocity of the TDGs in the 
M31 outskirts. 

We ado pt a baryonic -t o-tota l mass ratio of 20% 
as used bv lHammer et all l|2010h instead of 6% as in 
Yang & Hammer (2010), as well as in Sect. [4] Following 



Yang fc Hammerl |2010j), the dark matter scale length is 
tuned to fit the rotation curve. A total mass of 5.5 x 10 1 
M (3.3 x 10 11 M©) provides a R vtr = 70 kpc and c = 10 
(R vir = 60 kpc and c = 6) for M31 (the MW). This af- 
fects the lYang fe Hammerl (|20ld ) value for the M31 tan- 
gential velocity that is required for the LMC to be near 
M31 in the past, leading to (vx, vy) equals to (-135, 
10) km/s (fiw = — 71 /ias.yr -1 , /ijv = — 21 /ias.yr" 1 ) in- 
stead of (-105 , -7) km/s (jiw =_— 62 /ias.yr -1 , /ijv = 
—25 /ias.yr" -^. lYang fc Hammerl (|2010h show that such val- 
ues for the M31 tangential velocity are comparable to the 
radial velocity of M31 (129 km/s) and explain why they can 
be substantially larger than the tentativ e estimates made 
bv lvan der Marel fc Guhathakurtal l|200ct ). Our values (and 
Yang & Hammer values) ar e consistent within 1 . 7 (1. 1) a 
from the recent evaluation of lvan der Marel et alJ l|2012h us- 
ing very small proper motions of few M31 stellar fields based 
on a 5 year time basis. 

We first verify whether the adopted simple spherical 
potential for M31 in the analytic model is consistent with 
lHammer et al.l |2010l ) si mulations. After exami ning the dark 
matter distribution in the lHammer et all l)2010h simulations, 
we choose to pick up the tidal tail one Gyr after fusion, when 
the mass is relatively stable and spherically distributed. At 
this epoch, the MW is too far from M31 to affect the tidal 
tail. As this model does not allow us to follow the internal 
evolution of the tidal tail, we simplify the code by neglecting 
the force between the tidal tail particles. For each tidal tail 
particle the simulation is therefore a 3-body problem. In 
order to check whether the tid al tail has the same ev olution 
in the GADGET2 code used bv lHammer et alJ l|2010h that in 
the 3-body code used here, we ran two simulations with the 
same initial conditions. After 5 Gyr, the tidal tail positions 
and velocities are similar for the two simulations within 5%, 
which confirms that our choice for the initial conditions does 
not impact the results significantly. 

The different stages of the simulation are those de- 
scribed in Sect. 13.21 The tidal tail is selected from the GAD- 
GET simulation, one Gyr after fusion time, and inserted 
into the 3-body problem. This corresponds to 4 Gyr after 
the first passage. As the first passage is estimated to occur 
8.5±0.5 Gyr ago |Hammer et al.ll2010h . this corresponds to 
a look back time of 4.5±0.5 Gyr. 

5.3 Simulated VPOS properties 

The conditions for reproducing the VPOS using a merger 
model for M31 are considerably more constrained than in 



M31 

',-> 



MW trajectory 




Tidal Tail 
4.5 Gyr ago 



Figure 7. Schematic illustration of the VPOS formation as a re- 
sult of a tidal tail ejected by the M31 major merger. The tidal tail 
sweeps space for 4.5 Gyr followed by the LMC and finally passes 
through the MW trajectory. The green dots represent the initial 
and final position of the tidal tail, while white dots represents in- 
termediate steps. The M31 and MW disks are represented with a 
size of 100 kpc, as well as the observed VPOS (blue rectangular) 
and the y- and z-axis of the M31 rest-frame (green and red lines). 



Sect. [4] because of the required match between the tidal tail 
velocity and th at of the LMC. I n fac t, we realize that for 
most models in Hammer et al] (|2010h such a condition is 
verified as illustrated in Fig. [7] Besides the choice of the 
model reproducing the most accurately the properties of 
M31 (Yang et al., 2012, in prep.), the location of the tidal 
tail can be fine-tuned using a rotation by less than ± 20 
degre es along the M31 rotational axis (see lHammer et al.l 
2010). Such a fine-tuning preserves the modeling of the Gi- 
ant Stream and is required for matching the look-back time 
at which both the LMC and tidal tail are found at a common 
position in the M31 outskirts. Figure [7] shows the evolution 
of the tidal tail together with the LMC, and the MW tra- 
jectories relative to M31. It illustrates that for a given set of 
transversal velocities for the MW relative to M31, the LMC 
is found to travel between the two galaxies. Without more 
than fine tuning, to have its trajectory matched by the pre- 
dicted tidal tail formed after the first encounter between the 
two progenitors of M31. 

Let us now consider the first-passage tidal tail, which is 
the best candidate for inducing TDGs passing near the MW. 
Figure [5] shows the formation of clumps due to gravitational 
instabilities in the gaseous tidal tail at discrete locations. 
The five prominent over-densities have mass profiles consis- 
tent with that of dwarf galaxies with half-mass radii from 0.9 
to 4 kpc, and total masses from 1.7 to 8.6 x 10 s Mq (see Ta- 
ble |4}. Our 2.4 million particles simulation generates TDGs 
similar in location, e.g. the so-c alled "beads on a string " ap- 
pearance, to those found in the lWetzstein et al] (2007) EG2 
model with 6 million particles, but with larger masses. This 
is surprising, since the larger the resolution in the simula- 
tion, the larger t he mass of the formed T DGs (see Sect. 13.1] ). 
As suggested bv IWetzstein et aTl (|2007r ). this could be due 
to the initial gas fraction, whi ch is twice in ou r simu lation 
than the 30% adopted in the IWetzstein et all (|2007l ) EG2 
model and could result in the formation of more massive 
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Figure 8. Left panel: Gas density map of the first tidal tail, 9 Gyr after the first passage of the merger. Gas density is coded using 
blue-grey color, while stellar density is shown in orange. Five red circles indicate the location of the five most massive TDGs (mass larger 
than ~ 10 8 Mq), which are labeled TDGf to TDG5. Middle panel: Close-up of the TDG2. Right panel: Radial profile of the mass surface 
density of TDG2. The half mass radius ~ 1.5 kpc is indicated as a dashed line. 



TDGs. However, the simulation does not have a resolution 
large enough to sample the lowest mass dwarf galaxies (see 
Table Q. A higher resolution and a smaller softening lengt h 
is required to generate such TDGs (jBournaud et al.l 12008^1 . 
and one may expect a larger number of TDGs with masses 
similar to those of the MW companions. 

Given our absence of knowledge about the detailed in- 
ternal structure of the tidal tail, we have arbitrarily selected 
particles within a 40 kpc section of the tidal tail (300 parti- 
cles) at the region interacting with the MW (see red arrow 
in Fig. [HJ. In Fig. [9] for the sake of the visualization, we 
have randomly selected only 25 particles amongst the 300 
particles. They are all included into a thick plane, whose 
orientation is found to be similar to the VPOS within 13 
degrees, a value smaller than the observational uncertainty 
(see Table [2j . Finally, almost all their angular momentum 
directions are found to be close to the direction perpendic- 
ular to the observed VPOS, except for some of them which 
have an opposite direction as Sculptor (see Fig. I lOf) . 

However even if this tidal tail give some promising re- 
sults, it does not match the spatial distribution extension of 
the 11 classical dwarf galaxy (see Fig.[9l top panels). We no- 
tice that in the simulations reproducing M31 the first tidal 
tail may vary by factor 1 to 5 in width. This variation is due 
to the sensitivity to the initial conditions: the size of the 
gas disk in the progenitors and the orbital parameters. For 
example, orbital parameters associated to conditions close 
to the resonance (between the disk rotation and the orbital 
motion) at the first passage provide very narrow first tidal 
tails. Because being close to the resonance is not a require- 
ment to the M31 model (in fact such conditions are mostly 
required at the second passage), we have used another model 
of M31 to generate a wider tidal tail. As the clump locations 
within the tidal tail have no other constrains than the actual 
spatial distribution of the classical dwarf galaxies, we have 
used three clumps instead of one to mimic the extension of 
the tidal tail and the results are shown in Figure [9] (bottom 
panels) . 

We emphasize that this model is in sharp contrast with 



Table 4. Physical parameters of the five TDGs identified in 
the tidal tail generated during the 1st passage between the M31 
merger progenitors. Masses are given in unit of 10 8 Mq. From 
left to right: gaseous mass, stellar mass, total mass, gas fraction, 
and half mass radius (kpc). 



No. 


M gas 


Mstar 


^Total 


/gas 


rh 




(10 8 Mq) 


(10 8 Mq) 


(10 8 Mq) 




(kpc) a 


1 


1.67 


0.056 


1.72 


0.97 


4.3 


2 


6.20 


2.41 


8.62 


0.72 


1.5 


3 


3.49 


0.027 


3.76 


0.93 


1.1 


1 


4.11 


2.58 


6.70 


0.61 


0.75 


5 


4.09 


2.93 


7.02 


0.58 


0.85 



"r;, is the half mass radius of the TDG for the 2D projected mass 
profile. 

Sect. l4land |Pawlowski et al.l (|201ll ) studies that use all free 
parameters to optimize the tidal tail. In other words, the 
formation of the VPOS is predicted by the merger model of 
M31, as shown in Fig. [7] 



6 DISCUSSION 

6.1 A surprising series of "coincidences" 
supporting an M31 origin 

At first sight, a possible causal link between a possible 
merger in the history of M31 and the MW dwarf galaxies 
might appear at odds with many years of studies of the Lo- 
cal Group. To our knowledge such a proposition has never 
been made, probably because it appears far too unlikely for 
geometrical reasons: why the MW would be precisely at a 
rendez vous with a tidal tail formed during the first passage, 
8 to 9 Gyr ago, of a merger that occurred at more than 780 
kpc away? 

Let us consider three galaxies of similar masses, two of 
them being involved in a merger process. What is the chance 
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Figure 9. Face-on (left panels) and edge-on (right panels) views of the observed VPOS populated by the tidal tail particles (green 
diamonds). The 11 classical dwarf galaxies are represented by red points (dSph) and blue open squares (MCs). The top panel corresponds 
to the tidal tail with a small extension and the bottom one to the extended tidal tail. For both, there are 25 particles randomly selected 
into the clumps which interact with the MW (green diamonds). 



that a tidal tail produced during the event reaches the third 
galaxy? Assuming a merger in the history of M31, it should 
be rather small given the fraction of the volume swept by a 
tidal tail. The corresponding conic volume is shown in Fig. 
111! drawing a solid angle of 2000 square degrees, i.e, 5.5% of 
the 4-7T steradian sphere. In principle this argument should 
suffice to falsify or reject a possibility to associate M31 to 
the MW dwarf galaxies. In practice, the MW lies precisely 
within this small volume, because both the M31 disk is seen 



edge-on and the geometry of the Giant Stream (and of the 
Northern Loop) fixes the progenitors inclinations, and then 
the properties of the tidal tail induced after the first pas- 
sage. _We_jMve_coMidered a very large number of models 
from lHammer etafl (|2010h as well as more recent ones. In 
these models, r per i ranges from 24 to 30 kpc, the mass ratio 
from 2.5 to 3.5, the dark to baryonic matter ratio from 10 
to 25%, and a significant range of initial progenitor incli- 
nations are explored, which result in an uncertainty on the 
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Figure 10. Angular momentum directions of the selected particles for the small extension tidal tail (green asterisks) and the LMC (blue 
triangle). The dashed ellipse represents the uncertainty associated to the mean angular momentum of the 11 classical dwarf galaxies, 
while the solid ellipse represents the uncertainty associated to the direction perpendicular to the VPOS. 



angular location of the first tidal tail that is smaller than 
20 degrees. Of course we cannot claim having investigated 
all the numerous parameters of a merger model for M31, 
although the number of constraints just desc ribing the M31 
galax y and its outskirts are quite impressive l|Hammer et al.l 
|2010T ). We have also verified that f or all geometrical config- 
urations studied by I Barnes! (|2002T ) , the disk plane is close 
within less than 25 degrees to the orbital plane. This favor- 
able configuration applies well to the M31 disk (see Fig. llip . 
which is supported by an angular momentum that is found 
to be slightl y larger than the ave rage of similar local spiral 
galaxies (see lHammer et ai1l2007l ). 

Perhaps all of the above is just a coincidence and other 
properties would easily falsify a causal link between the M31 
merger and the MW dwarf galaxies. Fig. [7] shows that the 
tidal tail velocity matches well that of the LMC, in ampli- 
tude and orientatio n. This is ensured u sing the most recent 
results compiled bv lVieira et al.l (|2010l ) with quite large un- 
certainties (see their Fig. 10), while the amplitude of the 
LMC velocity would be too large using the lKallivavalil et al.l 
(2009) values. Such a match with the LMC velocity is an - 
other intriguing coincidence. While [Yang fe Hammerl (|2010l ) 
have shown that the LMC past trajectory may have passed 
near M31 for a precise range of M31 transverse velocity, why 
would the velocity of the modeled tidal tail follow such an 
amplitude and orientation? This coincidence naturally en- 
sures that our model can reproduce the VPOS orientation 
and dwarf angular momenta (see Fig. [9j . 

Possibly the most astonishing, fully unexpected coinci- 
dence is the fact that the part of the tidal tail shown in Fig. 
[8] which interacts with the MW is very close to the most 
massive simulated tidal dwarf, TDG2, whose mass almost 



reaches that of the LMC. This is very encouraging for pur- 
suing this study with higher resolution simulations. 



6.2 Could the MW dwarf galaxies be descendants 
of TDGs? 

6.2.1 Global properties 

First we acknowledge that a robust proof of a dominant 
dark matter component in MW dwarf galaxies would pro- 
vide a definitive falsification of the scenario proposed in this 
paper. It may come from the detection of dark matter anni- 
hilation in MW dwarf galaxies. In fact the non-detection of 
dark matter annihilation in MW dwarf galaxies could pro- 
vide the strongest up to date constraint on the nature of the 
dark matter particles, excluding for some annihilation chan- 
nels, thernT^Ld^xkjru^terfor masses below 30 GeV or so (see 
e.g.. lAckermann et al.|[2012T ). On the other hand, a tidal ori- 
gin for the MW dwarf galaxies would discard a significant 
dark matte r content in these g alaxies. It was conside red very 
early ( e.g., Lvn dcn-Bcll 197q) and was revived by iKroupal 
(|l997t) ; iKroupa et al.l (|2005l ) . Support s for this assumption 
are numerous and well summarised by IKroupa et alj (|2010l ) 
and lKroupal l|2012l ). Their old stellar population are perhaps 
inherited from the merger progenitors. Their long lifetime in 
the tidal tail could have let them in various environments 
providing a diversity in star for mation historie s that could 
explain their different behaviours. IKroupa et all (|2010T ) sum- 
marised that "the physics of TDG formation and evolution is 
sufficiently well understood to conclude that 1) once formed 
at a sufficient distance from the host, TDGs will take an 
extremely long time to dissolve, if at all; and 2) the TDGs 
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Figure 11. Geometrical view of the merger remnant. Top panel: 
The disk plane of M31 is represented using a picture of a spiral 
galaxy with its spin direction shown as a blue arrow, while the 
orbital plane of the progenitors as shown in green. The white cone 
represents the solid angle that could be swept by the tidal tail 
ejected during the merger. Bottom panel: Zoom on the solid angle 
swept by the tidal tail. The MW position and the plane of VPOS 
is shown in blue. The opening angle of the cone is constrained by 
the positions of TDG1 and TDG5 (see Fig. [8} and the modeling of 
M31 and the Giant Stream. It accounts for uncertainties related 
to the modeling (varying parameters, see the text) as well as for 
a rotation around the M31 disk by ±20 degrees along the M31 
disk and ±5 degrees in the plane of the tidal tail. 



formed will naturally lead to a population of ancient TDGs 
that resemble dSph satellites". 

However to be pl ausible our scenario (as well as that of 
IPawlowski et al]|201ll ) requires an alternative explanation of 
the lar ge mass-to-light rat ios of the MW dwarf galaxies. Re- 
cently Casas et al.l (|2012T ) have simulated dark-matter free 
satellites that after several orbits around the MW are out-of- 
equilibrium bodies with high apparent mass-to-light ratios. 
According to their claim, if such progenitors of MW satel- 
lites have reached the MW potential at the same time along 
a few specific orbits, they may reproduce most of their in- 
trinsic kinematics and surface brightness. Although similar 
simulati ons are beyond the scope of this paper, we noticed 
that the ICasas et al.l l|2012h progenitors are far more com- 
pact than the TDGs that are formed from our modelling. In 
other words, it is likely that realistic TDGs are more frag- 
ile and they may b e even more easily disrupted than in the 
ICasas et al.l l|2012h study. In a forthcoming paper we will 
study a more realistic TDG infalling into the MW potential 
and to verify whether or not such an event may reproduce 
the MW dwarf galaxies. Such a study is mandatory to ver- 
ify whether or not the large mass-to-light ratios of the MW 
dwarf galaxies falsify the M31 scenario. 



6.2.2 Could the LMC be descendant of a TDG? 

The above mentioned studies do not apply to the Magel- 
lanic Clouds and one may wonder whether or not the most 
massive Irr galaxies of the Local Group could also be TDGs. 
In the M31 scenario, both LMC and SMC have arrived for 
the first time in the MW halo, and as such they are not 
stripped of their gas, conversely to other MW dwarf galax- 
ies. The stellar mass of the LMC h as been properly esti- 
mated by Ivan der Marel et all (|2002h assuming an LMC V- 
band absolute magnitude of My = —18.5, an extinction of 
Ay=0A, and an exti n ction- corrected color (B — V)o = 0.43. 
Following iBell et al.l l^OOSt ), the stellar mass for the LMC 
is estimated to be 2.6 x lO 9 Af for a diet Salpeter IMF. 
Based on a full reconstru ction of the star formation history, 
lHarris fc Zaritskvl (2009) found a stellar mass ranging from 
1.7 to 3.5xl0 9 M© (see their Fig. 12) for a Salpeter IMF. 
To be compared wi th a diet Salpete r IMF value, one has to 
subtract 0.15 dex i|Bell et al.l [2003h to these values, reach- 
ing a range of 1.2 to 2.5 xl0 9 M Q . Other IMFs generally 
provides a lower ste llar mas s, by s ubtrac ting an add i tiona l 
0.15 dex for e.g. a iKroupal ([20021) or a iKennicuttl (Il983p 
IMF. According to lKim et al.1 l| 19981) . the LMC gas mass is 
0.5 x 10 9 Mq and thus the baryonic mass ranges from 1.7 to 
3.1 xl0 9 M Q , assuming a diet Salpeter IMF for the stellar 
mass estimate. These values are larger than those of Ta- 
ble |3] by a factor 2 to 3, which could be a serious difficulty 
for a tidal origin of the LMC. More detailed and resolved 
simulations of gas-rich tidal tail with more realistic physical 
conditions may he lp to further investig ate this issue. On the 
observational side. lKavirai et all (|201lh found several TDGs 
(approximately 15% of them) with a mass larger than 10 9 
Mq, after a wide TDG search in the SDSS data through the 
Galaxy Zoo project. The fact that these local mergers are 
certainly less gas-rich than the assumed ancient M31 merger 
progenitors suggests the latter may lead to TDGs as massive 
as the LMC. 

To investigate whether a tidal origin related to a merger 
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at the M31 location is realistic or not would require to com- 
pare it with the numerous results provided by studies on the 
Local Group and dwarf galaxies. It also requires to firmly as- 
sess the model properties and its predictions before any con- 
clusion. For example the recent determination of the Fe/H 
abundance of old stars in th e LMC may be not a falsifica- 
tion of this scenario (see e.g. lHaschke et al.ll2012h . The fact 
that the old globular clusters in M31 are relatively Fe/H 
richer can be explained by our model. The LMC is pre- 
dicted to be formed within the first passage tidal tail with 
stars stripped from the less massive interloper 8 to 9 Gyr 
ago, while old globular clusters could have been formed from 
material stripped from the most massive interloper, i.e., ex- 
pected to be metal richer. The search for further falsification 
tests has to continue. 



SIX) 




O l i i 

100 200 300 

Radius (kpc) 



6.2.3 MW dwarf velocities and angular momenta 

It is unclear yet whether all classical dwarf galaxies are part 
of the VPOS. Sagittarius is clearly off the VPOS in Fig. [3] 
Could it be due to the fact that it is strongly tidally dis- 
rupted by the MW potential? Figure [12] provides a crucial 
test of the M31 scenario by comparing the dwarf velocities 
to the escape velocities assuming two values for the MW 
mass. However this exercise is limited by the accuracy of the 
proper motion data: this excludes Draco, Sextans and LeoII 
for which the uncertainties are so large that their bound or 
unbound nature remains unknown. For example, the Sex- 
tans proper motion was not measured but derived, and the 
Draco proper motion was calculated by a ground-based tele- 
scope more than 10 years ago providing an extremely high 
value, V = 556 km/s. It results that Sag, UMi and Carina 
are certainly bound for any values of the MW mass, while all 
other dwarf galaxies are lying in a region where the measure- 
ment accuracy does not allow us to conclude. Perhaps the 
Sag, UMi and Carina properties are not consistent with the 
M31 merger scenario, and only much more accurate mea- 
surements of proper motion will allow to provide a similar 
test for the other dwarf galaxies. Using proper motions to 
distinguish between the various scenarios for the formation 
of the dwarf galaxies is ver y tempting but ofte n limited by 
the data accuracy (see e.g., I Angus et al.ll201lh . and by the 
poor knowledge about the true MW total mass. Moreover 
this also depends on the nature of their progenitors. If as- 
sumed to be of tidal origin from a gas-rich merger, they 
should be gas-rich as shown in Table [4] prior they interact 
with the MW potential. 

Possibly, when entering the MW potential, these gas- 
rich TDGs are entirely stripped of their gas through 
tidal shocks and r am pressure stripping as described by 
Maver et all j2007Ti . More recently, many authors (see e.g., 
Maverl 2010l . 1201 ll . and references therein) have shown that 
tidal stirring is an efficient mechanism to remove the gas 
from gas-rich, dark matter dominated, dwarf irregular galax- 
ies interacting with the potential of a massive galaxy such 
as the MW. This leads to their transformation into dwarf 
spheroidals, which turns rotationally supported stellar sys- 
tems into pressure supported ones. The interplay between 
dynamical friction and tidal mass evaporation not only de- 
creases strongly their mass by factors reaching 10, it af- 
fec ts also their angula r momentum and velocity as shown 
bv lTaffoni et~afl (|2003h . This might lead to a ~ 20% decay 



Figure 12. MW escape velocity as a function of radius for two 
baryonic-to-total mass ratios: 6% (magenta dashed line) and 20% 
(blue line). The baryonic mass remains equal to 6.6xl0 10 M 
in both cases. The red lines represent the uncertainty associated 
with the velocities of the classical dwarf galaxies. 

in velocity during 2 Gyr (see their Fig. 6), e.g., providing 
initial velocities of UMi and Carina consistent with those 
of other dwarf galaxies. Our initial guess is that tidal stir- 
ring should be more efficient for progenitors that are not 
dark matter dominated such as TDGs. These are however 
pure conjectures that underlines the severe limitations of our 
modeling for which we have represented tidal tail particles 
as point masses. For reaching a conclusive answer about the 
consistency of the M31 merger scenario with the dwarf prop- 
erties would require a far more complete modeling, perhaps 
taking into account the interaction between the whole tidal 
tail and its associated TDGs with the MW potential. Up 
to our knowledge, such a simulation have never been done 
by assuming a dark-matter free irregular galaxy such as a 
realistic TDG falling in the potential of a massive galaxy. 
We are encouraged to perform this in a future study also 
because we notice that the most inconsistent galaxy veloc- 
ities (those of Sag, UMi and Carina) have also the angular 
momenta the most shifted from the averaged dwarf galaxy 
angular momentum (see Fig. [3} . 



7 CONCLUSION: A CAUSAL LINK BETWEEN 
THE M31 MAJOR MERGER AND THE 
VPOS? 

It is widely accepted that the class ical bulge of M31 is 
due to an ancient major m erger l|van den Berghl 120051 : 
iKormendv fc Kennicuttl |2004) . Quoting John Kormendy 
(2010, private communication): "don't we already know with 
some confidence that M31 was formed in at least one (in 
fact: very probably exactly one) major merger?" Perhaps 
its consequences have not been fully investigated, whereas 
it occurred in the galaxy that represents almost two third of 
the Local Group baryonic mass. The immediate question is: 
where are the relics of the most energetic event that occurred 
in the whole Local Group history? One may investigate the 
most anomalous and difficult to explain features discovered 
in the Local Grou p. Among them are the M31 haunted halo 
l|lbata et alj2007l ). the Giant Stream (jlbata et alj200ll ). the 
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VPOS and also the u nusual proximity of t wo massive Irr 
near the MW (see e.g.. Ivan den Berghll2010l ). 

In principle a ~ 1 Gyr old minor me rger may reproduce 
the Giant Stream (e.g.lFardal ct al. 2008) although this may 
be problematic (see e.g. Font et al]|2008l ) bec ause stars in 
the G i ant Stream have ages from 5.5 to 13 Gyr (|Brown et al.l 
l2007l ). lHammer et all (|2010h have been able to reproduce the 
Giant Stream as well, assuming a major merger and using 
the stellar ages in the M31 outskirts as clocks to date a ma- 
jor merger, providing a first passage and a fusion occurring 
8.75 ± 0.35 Gyr and 5.5 ± 0.5 Gyr ago, respectively. Having 
M31 involved in a merger at these epochs may be quite com- 
mon, as demonstrated by detailed studies of intermediate 
redshift galaxies from t he IMAGES sample l|Hammer et al.l 
120091 ; iPuech et al.1120121 ). This model is quite challenging in 
reproducing many other properties of the M31 galaxy (disk, 
bulge, thick disk and 10 kpc ring) together with the Giant 
Stream. A noticeable feature of gas-rich major mergers is 
that they produce long-lived tidal tails that can reproduce 
many struct ures discovered in the halo of giant spiral galax- 
ies (see e.g. JWang et al.ll2012l ). 

This paper is the first of a series aiming at investigate 
consequences of such a major merger occurring in the Local 
Group. The main weakness of the present modeling is due 
to the oversimplification in considering tidal tail particles as 
point masses. To reproduce the properties of the MW dwarf 
galaxies requires a full modeling of extended, gas-rich TDGs 
entering into the potential of the MW. It is highly desirable 
to verify whether or not the MW dwarf galaxies and the 
VPOS can be fully modeled by the interaction between a 
gas-rich tidal tail, its associated TDGs, and the MW. Such 
a study will aim at verifying whether or not the intrinsic 
properties of dwarf galaxies may be reproduced, including 
their apparent large mass-to-light ratios. 

Our modeling of t he M31 major merger is getting more 
and more mature since lHammer et al.l l|2010l ) and it is able to 
provide some predictions that can be tested for falsification 
purposes. The modeling of an anisotropic feature like the 
Giant Stream, and further on, of the Northern loop (Yang 
et al. 2012, in prep.) almost fix the location and velocity of 
the tidal tails associated to the merger. Here we have con- 
centrated our efforts to examine the largest tidal tail, which 
is formed during the first passage, stripping the minor en- 
counter stars, 8 to 9 Gyr ago. This tidal tail is bringing the 
largest amount of baryonic material expelled by this event 
towards the Local Group. It results in a series of "coinci- 
dences" that considerably strengthens our preliminary con- 
jecture. The MW is found precisely at the meeting point 
with the tidal tail, at present time. Moreover, the velocity 
of the tidal tail matches well that of the LMC providing a 
good fit of the VPOS, including its spatial orientation and 
the angular momentum, even to some details such as the 
inverted angular momentum of Sculptor. 

Having this in mind we conclude that a link between 
the VPOS and the major merger at the M31 location is 
plausible. Compared to the MW merger or encounter hy- 
pothesis of iPawlowski et ail ([20 11), it has the advantage of 
resulting from a testable prediction that could be falsified 
or supported by many observations of the Local Group. Be- 
sides this, the present modeling is made within the ACDM 
scenario and there is no peculiar need to account for other 
physics to describe the above noticed unusual properties of 



the Local Group, i.e. the VPOS, the presence of two mas- 
sive Irr near the MW and most features found in the M31 
halo. In fact, iKnebe et all (|201lT ) have been the first to show 
that within ACDM, there should be renegade sub-halo, i.e. 
sub-halo of M31 that could have migrated towards the MW. 
However they could not explain neither the VPOS nor the 
LMC and its large velocity, conversely to our prediction that 
lead to an external origin for most MW dwarf galaxies. If 
the latter is true, it strengthens the tension between the pre- 
dicted number of dark matter sub-haloes and observations, 
i.e. the so-called missing satellite problem, which might be 
then reconsidered as an excess of sub-haloes in the ACDM 

paradigm. 

Note added into proofs : IPawlowski et a l. (2012b) re- 
cently shows that filamentary infall is excluded to explain 
the VPOS. 



ACKNOWLEDGMENTS 

This research has made use of the NASA/IPAC Extragalac- 
tic Database (NED) which is operated by the Jet Propul- 
sion Laboratory, California Institute of Technology, under 
contract with the National Aeronautics and Space Admin- 
istration. We are grateful to Jianling Wang for generating 
a part of the models of M31 that have been used in this 
study. Finally, to represent the Local Group in 3D, three- 
dimensional visualization was conducted with the S2PLOT 
programming library. We warmly thanks the referee whose 
remarks have significantly improved the manuscript includ- 
ing by identifying the important tests that should be made 
to put the M31 scenario in a cosmological context. 



APPENDIX A: STATISTICAL METHOD TO 
TEST THE ISOTROPY 



The statistical test proposed by I Met z et al.l |2007al ) assesses 
whether the MW dwarf galaxy positions or velocities de- 
rive from an isotropic distribution. 50 000 random samples 
are generated using Monte-Carlo simulations. Each sam- 
ple is isotropic and follows a radial density distribution 
p{r) oc 1/r 2 , consistent wit h the radial distribu tion of the 
MW dwarf galaxy positions |Kroupa et al.| [2005) and veloc- 
ities. A plane is fitted for each sample, and the associated 
uncertainty is estimated using bootstrap resampling as fol- 
lows. Each sample is resampled 1 000 times and the perpen- 
dicular vector for each of the corresponding fitted plane is 
stored in a matrix M. The 3D symmetric matrix S = M T M 
is computed and its three eigen values (n < T2 < T3) cal- 
culated. To characterize the 1 000 vector distribution, i.e., 
to characterize the plane fit uncertainty, two numbers are 
defined using the eigen values: 

= log (t 3 /t 2 ) 
log (t 2 /ti) 

C = log (r 3 /ri) 

7 describes how clustered the vector distribution is, while 
£ indicates how the vector distribution focuses towards one 
direction. If 7 and £ are small (~ 1) the vector distribution 
is close to isotropic (the plane fit uncertainty is thus large), 
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whereas if 7 and C are large (> 2) the vector distribution 
clearly defines one direction (the plane fit is then well de- 
fined). The resulting statistics of 7 and C is plotted in Fig. [3] 
For the f classical dwarf galaxy position and velocity set, 7 
and C are also derived. The distribution function D(7o, Co), 
defined by the probability that 7 > 70 and C > Co, can be 
used to estimate the probability for a sample to derive from 
an isotropic law. 
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